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Abstract 

We investigate spherically symmetric continuously self-similar (CSS) solutions in the SU{2) 
sigma model coupled to gravity. Using mixed numerical and analytical methods, we provide 
evidence for the existence (for small coupling) of a countable family of regular CSS solutions. 
This fact is argued to have important implications for the ongoing studies of black hole 
formation in the model. 

1 Introduction 

Self-similar solutions of Einstein's equations have been extensively studied in general relativity 
(see H] for the recent review and vast bibliography of this subject). They are interesting for 
several reasons. First, under the assumption of self-similarity Einstein's equations simplify con- 
siderably and sometimes (in spherical symmetry, for instance) reduce to a system of ordinary 
differential equations. This enables one to study them using dynamical systems methods and 
in some cases it is even possible to find explicit solutions. Second, self-similar solutions exhibit 
a well-defined singular behaviour at the blow-up. This makes them relevant in the studies of 
singularity formation; in fact, most of the examples of naked singularities in the literature involve 
self-similar solutions. Third, in some situations solutions of the Cauchy problem starting from 
generic initial data evolve asymptotically to self-similar form, in other words, self-similar solu- 
tions can play the role of attractors. Finally, the recent surge of interest in self-similar solutions 
is due to the fact that they appear as critical solutions at the threshold for black hole formation 
in gravitational collapse. 

Such reasons motivated also this paper which, in a sense, is an extension of Q where self- 
similar wave maps from Minkowski spacetime into the 3-sphere were studied. It was proved 
in that there exists a countable family of regular (by regularity we mean analyticity below 
the Cauchy horizon of the singularity) CSS solutions, labeled by a nonnegative integer n (a nodal 
number). It was also shown that the nth solution has exactly n unstable modes. The role of 
these solutions in the dynamical evolution was investigated in |^ where it was shown that: (i) 
the n = solution determines a universal asymptotics of singularity formation, (ii) the n = 1 
solution appears as a critical solution at the threshold for singularity formation (see also 0). 
Here, we generalize by turning on gravity, that is, we consider wave maps with the 3-sphere as 
the target and a domain manifold which is not Minkowski but a spacetime satisfying the Einstein 
equations with the SU{2) sigma field as the source. The coupling of gravity is parametrized by 
a dimensionless coupling constant a (so that a = corresponds to gravity being turned off). 
Since the coupling does not break scale invariance, it is natural to ask if the a — Q CSS solutions 
constructed in l2|, persist for nonzero a. A positive answer to this question is the main result of 



our paper. More precisely, we show that the nth CSS solution which is analytic within the past 
light cone of the singularity (in what follows we will refer to the past and future light cones of 
the singularity as to the past and future self-similarity horizons (SSH)), persists up to a = 1/2. 
However, this solution is regular (that is, analytic up to the future SSH) only for a < an, where 
{an} is an increasing sequence bounded above by 1/2. For a„ < a < 1/2, the past SSH of the 
nth solution is surrounded by a spacelike apparent horizon beyond which the solution cannot 
be smoothly continued. Although we have not studied linear stability, it seems plausible that 
the stability properties of solutions for small a are the same as for a = 0. Assuming this, the 
n — solution (for a < ag ~ 0.0688) is a generic naked singularity, while the n = 1 solution 
(for a < ai ~ 0.1518) is a codimension-one naked singularity and the candidate for a critical 
solution. 

The CSS solutions constructed here should be relevant for the character of type II critical 
collapse in the model which is currently being investigated by Husa, Lechner, Piirrer, Thorn- 
burg, and Aichelburg (HLPTA). Their first report of progress Q shows discretely self-similar 
(DSS) behaviour at the threshold for black hole formation for a > 0.18 and indicates that DSS 
disappears for small a. It is natural to expect that for a < ai the n — 1 CSS solution takes 
over as a critical solution. Preliminary results of HLPTA seem to confirm this expectation [^. 
A mechanism of the DSS/CSS changeover is not yet understood, however. We anticipate that 
the standard type II CSS critical behaviour will be observed only in the range ao < a < ai. For 
a < ao we expect qualitatively the same behaviour as in the flat spacetime, where the n = 1 
solution sits at the threshold for the n = naked singularity formation. 

The rest of the paper is organized as follows. In section 2 we define the model and derive the 
field equations. In section 3 we construct a countable family of CSS solutions which are analytic 
between the center and the past SSH. An extension of these solutions beyond the past SSH is 
determined in section 4. In section 5 we speculate about the role of regular CSS solutions in the 
dynamical evolution. 



2 Field Equations 

Let X : M ^ N he a. map from a spacetime {M,gab) into a Ricmannian manifold {N,Gab)- 
Wave maps coupled to gravity are defined as extrcma of the action 



M 



with the Lagrangian density 



LwM - -^g'^'daX^dhX^'GAB- (2) 



Here G is Newton's constant and is the wave map coupling constant. Remarkably, in the 
"physical" case of 3 -I- 1 dimensional spacetime (to which our paper is confined), the coupling 
of gravity does not break scale invariance but introduces the dimensionless coupling constant 
a = AnGf^. The field equations derived from (1) are the wave map equation 

DgX^ + TiciX)daX^di,X^g-'' = 0, (3) 

where r^(j{X) are the Christoffel symbols of the target metric Gab and Dg is the d'Alembertian 
associated with the metric gab, and the Einstein equations Rab—\gabR — SnGTab with the stress- 
energy tensor 



Tab = fl [daX^dbX^ - ^-gab{g"'d,X''ddX'')^ G AB ■ 



(4) 
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The target manifold is taken as the three-sphere with the standard metric in polar coordinates 

= {F, e, $) 

GABdX^dX^ = dF^ + sin^F {dQ^ + sin^e (5) 

For the domain manifold we assume spherical symmetry and use Schwarzschild coordinates 

Qabdx^dx^ = -e-'^^A d^ + A-^dr'^ + r^{de'^ + sin'^e dcj)^), (6) 

where 5 and A are functions of {t,r). Next, we assume that the wave maps are corotational, 
that is 

F = F{t,r), Q = e, ^ = ^. (7) 
Eq.(3) reduces then to the single semilinear wave equation 

where 



□g = -e'dt{e'A-'dt) + ^dr{r^e-^A dr), (9) 



and the Einstein equations are 



dtA = -2arA{dtF){drF), (10) 
dr5 = -ar{{drFf + A-^e^\dtFf), (11) 

drA = ^-ar(A{drFf+A-^e'\dtFf + 2'^). (12) 



These equations are invariant under dilations (t, r) (At, Xr) so it is natural to look for con- 
tinuously self-similar (CSS) solutions, that is solutions which are left invariant by the action of 
the homothetic Killing vector K = tdt + rdr- To study such solutions it is convenient to use 
similarity variables p = r/{—t) and r = — ln(— t). Then K = —dr, so CSS solutions do not 
depend on r. Assuming this and using an auxiliary function z = e^p/A, we reduce Eqs.(8-12) 
to the system of ordinary differential equations (where prime is d/dp) 

F" + V-„(l+.>F"-3i=g^=0, (13) 



A' = -2apy4F'^ (14) 
pz' = z{l + a{l-z^)p^F'^), (15) 
pA' = 1- A- a (^p''A{l + z^)F'^ + 2 sin'' Fy (16) 

The combination of (14) and (16) yields the constraint 

1- A- 2a sin' F + aAp^F'^ {I- z^) = 0. (17) 

In what follows, we will be solving Eqs. (13-15). The constraint (17), which is the first integral 
of these equations, could be used to eliminate A from (13) but we find more convenient not to 
do so. The system (13-15) has a fixed singularity at the center p = and moving singularities 
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at points where z(p) ~ ±1 and/or A(^p) = 0. In terms of the similarity coordinate p, the metric 
(6) takes the form 

ds^ = A-^(\ - 2^2) p^dt^ + 2A-^tp dtdp + A-^t^dp^ + t'^p'^ide^ + sin^e dcfP), (18) 

hence the hypersurfaces z = ±1 are null (provided that A> 0). After we will refer to these 
characteristic hypersurfaces as to the self-similarity horizons (SSH). The first po where z(po) = 1 
is the locus of the past light cone of the singularity at the origin {t = {),r — 0) and we will call it 
the past SSH. We use the remaining coordinate freedom to locate it at po = li that is z(l) = 1. 
In the next section we will analyze the system (13-15) below the past SSH, that is for < p < 1. 

3 From the center to the past self-similarity horizon 

Expanding about p = 1, we get a one-parameter family of solutions of Eqs. (13-15) 

F(p)^ |+6(p-l), z(p)^l + (p-l), A{p)^l-2a-2a{l-2a)b'{p-l). (19) 

We shoot these initial data towards p = and adjust the parameter h so that the solution 
satisfies the regularity condition at the center, 

F(0) = 0, z(0)=0, ^(0) = 1. (20) 

We find numerically that for each a < 1/2, there is an infinite sequence {bn}, n = 0, 1, . . . , such 
that the corresponding solutions {Fn, An, Zn) satisfy the boundary conditions (19) and (20). The 
index n denotes the number of solutions of the equation Fn{p) = 7r/2 on the interval < p < 1. 
In the limit a — > 0, the nth solution tends uniformly to the corresponding a = solution. In 
the limit a 1/2, the parameters 6„(a) diverge to infinity indicating that regular solutions 
disappear for a > 1/2. Fig. 1 shows how the shooting parameters 5„ depend on a for n = and 
n = 1. In Fig. 2 we plot the n = and n — 1 solutions for a — 0.1. 

To summarize, we claim that for each a < 1/2 there exists a countable family of CSS solutions 
which are analytic at the center and at the past SSH. The structure of these solutions is basically 
the same as for a = 0. Qualitatively new things happen beyond the past SSH and now we turn 
our attention to this much more interesting region. 

4 Beyond the past self-similarity horizon 

For a = it was easy to show that all solutions starting at p = 1 with initial values (19) 
remain regular for all p > 1. As we shall see below this is still true for a > provided that 
b is sufficiently small (note, however, that the metric function A{p) is monotone decreasing, so 
A{oo) < 1, which means that the geometry is not asymptotically flat but conical). However, if b 
is large, the function z{p) is not monotone increasing for all p > 1 and a "sonic point" develops 
for some finite ps > 1, that is hmp^p^ z{p) — 1. We stress that this "sonic point" is not a SSH 
but an apparent horizon because limp^p^ A{p) = as follows easily from (14) and (15). These 
two kinds of behaviour are illustrated in Fig. 2. 

Note that p — oo corresponds to the hypersurface (t = 0, r > 0) so in order to analyze the 
global behaviour of solutions (for t > 0) we need to go "beyond p — oo". To this end we define 
a new coordinate x by 

d d , ^ , , 

-=p.-, .(p^l)^O. (21) 
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We also define an auxiliary function w{x) = 1/ z{p). In these new variables, the past SSH where 
ui = 1 is at a; = 0, while the future SSH (which nota bene is a Cauchy horizon) is located at 
some > where w{xa) = — 1- 

In terms of x and w, Eqs. (13-15) become autonomous (where now prime is d/dx) 

w' = -l + a(l-w2)i^'^ (22) 
A' = -2aAwF'^, (23) 

iAF'Y = (24) 

The constraint (17) becomes 

l-A-2asm'^F + aAF'^{w^ -1)^0. (25) 

From (19) the initial conditions at a; = are 

P{x)^^+bx, w{x)^l-x, A{x) ^ l-2a-2a{l-2a)b^x. (26) 

We already know from the previous section that for each a < 1/2 there is an infinite sequence 
{&„} determining solutions which are regular for all x < (note that p — corresponds to 
X = — oo). Although these solutions are completely fixed at the past SSH by the requirement of 
analyticity at the center, in order to understand their behaviour for x > it is helpful to drop 
this requirement temporarily and consider solutions that start from generic initial conditions 
(26). Of course, such solutions in general arc not analytic at the center. 

We first shall show that if b is sufficiently small, then the solution remains regular up to the 
future SSH. To see this, let / = (F - 7r/2)/fe. Then, in the hmit 6^0, Eqs. (22-24) reduce to 

w' = -l + a{l-w^)b^f'^^-l, (27) 
A' = -~2aAwb^f'^ ^0, (28) 

. ^^^^^ 

b(w — 1) w — 1 

with the initial conditions 

/(0) = 0, /'(0) = 1, w{0) = l, A{0) = l-2a. (30) 

The limiting equations (27) and (28) are solved by w = 1 — x and A = 1 — 2a. Substituting 
these solutions into (29) we get the equation 

(l-2.)r + -^^0, (31) 

whose solution is given by the hypergeometric function /(x) — X2Fi[{l—i^/7)/2, {l+i^/7)/2,2,x/2]. 
Since solutions of Eqs. (22-24) are continuous in b and x, by uniform continuity on compact in- 
tervals, the solutions with sufficiently small b will tend to w = — 1. 

Asymptotic analysis at the future SSH, xa, yields the following leading order behaviour 
(where y — xa ^ x) 

w^-l + y, A^ Ao-2aAoc^yh-i\y), F ^ F^ - cy\n{y), (32) 

where Aq ~ \ — 2asin^-Fb7 c = — sin(2Fo)/2ylo, and Fq is a free parameter. Thus, a solution 
that is analytic at the past SSH, generically will be only at the future SSIl|^. Nevertheless, 

^If -Fo is an integer multiple of 7r/2, then the singular log terms in (32) are absent and the solution is smooth 
at the future SSH. One can show that this happens for discrete values of b, which means that there exist solutions 
which are analytic both at the past and the future SSH. Of course, for generic a these solutions are singular at 
the center. 
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we checked that with the asymptotic behaviour (32) all curvature invariants remain finite as 
X — > xa ■ We interpret this somewhat surprising fact as an indication that the singularity at the 
origin is nakedQ. 

Next, we shall show that solutions with large b develop a "sonic point" for some finite xb > 0, 
that is, Imij:—,^^ w{x) = 1 and limx^xB ^i^) — 0- This time, we define the variables 

e = h{0 = b^{l-w{x)), s{0=b{F{x)~^). (33) 

Then, in the limit b oo, Eqs. (22-24) reduce to (where now prime is d/d£,) 

h' = 1 - a/i (^2 - ^ j s'^ ^ 1 - 2a/^s'^ (34) 

A' = ~2aA(^~ ^s'^ -2aAs'^, (35) 

with the initial conditions 

h{Q) = 0, ^(0) = 1 - 2a, s(0) ^ 0, s'(0) = 1. (37) 

It follows from (36) and (37) that, in the limit b oo. As' = 1 — 2a. Plugging this into (34) 
and (35), and using (37), we get the limiting solution 

AiO = (l-2a)Vl-4< and /i(0 = 7^ Vl " 4< (1 - ^1 - 40- (38) 

2a 

This solution develops a "sonic point" at ^ = l/4a. Again, by uniform continuity on compact 
intervals, we conclude that for solutions of Eqs. (22-24) with large 6 (and nonzero a), the function 
w{x) attains a minimum and then tends to 1 at some x —^ xb, while the function A(x) drops to 
zero at a;B. The leading order behaviour at the apparent horizon at x_b is 

w ^ 1 ~ a^/xB ~ X, A ^ d^/xB — X, F ^ Fb H ^^/xb — x, (39) 

where the positive parameters Fb ,a,d are constrained by the relationship 1 — 2ck sin^ Fb — ad/2 = 
0. Substituting this asymptotic behaviour into the metric (18), one can readily verify that the 
hypersurface x = xb is spacelike. 

Thus, solutions which start from the past SSH with initial conditions (26), tend in finite 
"time" to w = —1 if 6 is small, or to w = +1 if b is large. In what follows, we will refer to 
these two behaviours as to type A and type B solutions, respectively. The solutions of type A 
obviously form an open set and we believe (but have no proof yet) that the same is true for the 
solutions of type B. If so, there must exist solutions which are not of type A or B, call them type 
C solutions. By definition, type C solutions remain in the strip w € (—1, 1) which implies that 
they exist for all x (because one can easily show that solutions can go singular only if w — > ±1)- 

Our numerical analysis of the structure of types A, B, and C can be summarized as follow^ 
For a given a, there exists a critical value b*{a) such that solutions with b < b* are of type A, 

similar behaviour for CSS solutions of the Einstein-axion-dilaton equations was observed by Eardley, 
Hirschmann, and Horne |8|. 

^To simplify the exposition, we are cheating at this point. We have evidence that for large values of a (~ 0.42) 
the transition between types A and B occurs not at a single point but in a narrow interval (b^j^, femax)- We 
suspect that this interval, which can be thought of as the boundary between the basins of attractions of type A 
and B behaviours, has a fractal structure. This intriguing fact, suggesting perhaps a chaotic behaviour, is under 
investigation and will be described elsewhere. 
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solutions with b > b* are of type B, and the solution with 6 = 5* is of type C. Thus, we have 
a bistable behaviour with two generic final states A and B, and the separatrix C. Although the 
precisely critical initial condition b = b* cannot be prepared numerically, in order to figure out 
a behaviour of the type C solution it is sufficient (due to continuous dependence of solutions 
on initial conditions) to determine the flow of nearly critical initial conditions b — b* ± e. As 
is shown in Fig. 3, before the marginally critical type A and B solutions tend to w — ±1, they 
exhibit a long transient almost periodic behaviour. This leads us to conjecture that the type C 
solution is asymptotically periodic. The periodic intermediate attractor is clearly seen on the 
(w. A) plane. 

The complete discussion of periodic solutions would take us too far afield, so we postpone it 
to the subsequent paper. Here, let us only point out that the behaviour of CSS solutions outside 
the past SSH resembles some aspects of the DSS type II critical behaviour at the threshold for 
black hole formation. Actually, it can be interpreted as a poor ODE version of this phenomenon. 
The type A and B solutions are the analogues of dispersive and black-hole solutions, respectively. 
The type C solution is the analogue of the critical evolution which tends asymptotically to the 
critical periodic solution (which is the analogue of the choptuon). We even have the analogue of 
the black-hole-mass scaling. Namely, let A be the unstable eigenvalue (that is Re{X) > 0) around 
the periodic solution. Then the distance of a nearly critical flow (corresponding to b — b* + e) 
from the separatrix C is proportional to eexp(i?e(A)x), which means that the length of the 
transient periodic phase is proportional to — ^^\^\- locus of the future SSH xa (for the 
marginally critical type A solutions) and the locus of an apparent horizon xb (for the marginally 
critical type B solutions) scale in the same manner (see Fig. 4). This leading order scaling law 
is decorated in the next order by periodic wiggles with period Re{X)T. The origin of these 
wiggles is basically the same as in the case of black-hole mass scaling in the type II DSS critical 
collapse [0 . In passing, we remark that the periodic solutions are nonperturbative in a ~ in the 
limit a 0, the period and the amplitude of oscillations go to zero. 

Now, we return to the countable family of CSS solutions constructed in section 3 and finally 
determine their extension beyond the past SSH. As follows from the above discussion, for a 
given a, the behaviour of the nth solution outside the past SSH depends on whether is 
smaller or larger than b*: the solutions with 6„ < b* are of type A (and therefore regular in our 
terminology), while the solutions with 6„ > b* are of type B. The function b*{a) is monotone 
decreasing and b*{a) — > oo as a — *■ 0, while the functions fo„(a) are monotone increasing, hence 
we obtain an infinite sequence {a„} defined by the intersections &*(«„) — bn{an) (see Fig. 1). 
Since &„+i(a) < 6„(a) and bn{a) ^ oo as a — > 1/2, the sequence {a„} is increasing and bounded 
above by 1/2. Therefore, we conclude that the nth CSS solution from section 3 is regular iff 
a < an- The first two values of a„ are ao ~ 0.0688 and ai « 0.1588. 

5 Final remarks 

The role of the CSS solutions considered in this paper in the dynamical evolution depends 
crucially on their stability with respect to small perturbations. Although we have not performed 
the stability analysis, it seems plausible that for sufficiently small a the stability properties are 
the same as for a = 0, that is the nth solution has exactly n unstable modes. Speculations that 
follow are based on this assumption. We expect that for a < ao the situation is qualitatively the 
same as in the flat spacetinie, that is the n = solution determines the universal asymptotics of 
a naked singularity formation and the n = 1 solution sits at the threshold of this process. For 
a > ao the singularity of the n = solution is covered by a horizon so generic naked singularities 
presumably disappear - nota bene this could be interpreted as a gravitational desingularization 
phenomenon. In the range a S {ao,ai) we expect the type II threshold behaviour with the 
n = 1 CSS critical solution. What happens for a > ai « 0.1518? Since the n = 1 solution is 
not regular, it cannot be a bona fide critical solution (it can still appear as a local intermediate 
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attractor, though). HLPTA showed in Q that for large values of a the critical behaviour is DSS 
and the echoing period A increases rapidly as a tends to 0.18 from above (see Fig. 4 in |^). 
Moreover, they noticed small deviations from exact DSS at this smallest a value for which they 
observed DSS critical collapse. These facts led them to conjecture that a critical solution ceases 
to be DSS for still smaller values of a. This conjecture fits nicely with our conjecture that a 
critical solution ceases to be CSS for a > ai. Taken together, these conjectures suggest that the 
transition from CSS to DSS critical behaviour is not sharp but rather has a crossover character. 
HLPTA are currently running the computer simulations of the critical behaviour for the coupling 
constants a < 0.18, so the above speculations should be verified (or falsified) soon. Hopefully, the 
numerical phenomenology will also help understand a yet unknown mechanism of the CSS/DSS 
changeover. 
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Figure 1: The shooting parameters for n = and n = 1. The critical curve b*{a) 

separates solutions that have an apparent horizon beyond the past SSH from those that do not. 
The intersections of the curve b*{a) with the curves 60(0;) and 61(a) determine the critical values 
of the coupling constant ao ~ 0.0688 and ai « 0.1518. 
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Figure 2: The n = and n = 1 CSS solutions foid; = 0.1. Since ao < 0.1 < ai , the n = 1 solution 
exists for all p > 1, whereas the n = solution develops an apparent horizon at w 3.3683. 




Figure 3: The metric functions w and A for tMo nearly critical solutions for a — 0.2: type A 
solution with b = 0.36425022396604 (dashed line) and type B solution with b = 0.36425022396605 
(solid line). For intermediate values of x the solutions are almost periodic with the period 
T Ri 1.418. The intermediate attractor is seen on the {w, A) plane as an unstable limit cycle. 
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Figure 4: The locus of the future SSH. xa. against the logarithmic distance from the critical 
value, ln(&* — 6), for a = 0.2 is plotted as the solid line. The dashed line (shifted down for better 
comparison) shows the least-square fit of the leading order behaviour xa = — ^^i^* ~ b) with 
Re{X) « 2.029. The wiggles superimposed on this straight line are the imprint of periodicity of 
the intermediate attractor. The numerically computed period of the wiggles is approximately 
equal to 2.887 which agrees with the theoretically predicted value i?e(A) T. 
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